Intrinsic avalanches and collective phenomena in a Mn(II)-free radical ferrimagnetic 

chain 
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Magnetic hysteresis loops below 300 mK on single crystals of the Mn(II) - nitronyl nitroxide free 
radical chain (Mn(hfac)2(-R)-3MLNN) present abrupt reversals of the magnetization, or avalanches. 
We show that, below 200 mK, the avalanches occur at a constant field, independent of the sample 
and so propose that this avalanche field is an intrinsic property. We compare this field to the 
energy barrier existing in the sample and conclude that the avalanches are provoked by multiple 
nucleation of domain- walls along the chains. The different avalanche field observed in the zero field 
cooled magnetization curves suggests that the avalanche mechanisms are related to the competition 
between ferromagnetic and antiferromagnetic orders in this compound. 

PACS numbers: 75.50.Xx, 75.10.Pq, 75.60.-d, 75.60.Jk 
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I. INTRODUCTION 

Magnetic avalanches, that is to say abrupt magnetiza- 
tion reversals, often appear in low temperature hysteresis 
loops*!* 2 -! 3 -. In some cases, they are just giant Barkhausen 
jumps and so give information about defaults in the sys- 
tem. An interesting situation arises when avalanches are 
intrinsic to the system. Then they can be a powerful 
tool for understanding the underlying interactions, and 
the mechanisms of nucleation and domain-wall propaga- 
tion. 

Much interest has been devoted to the study of spin dy- 
namics in molecular nano-cluster zero-dimensional (OD) 
systems and one-dimensional (ID) or quasi-lD systems. 
In the latter, because of the strong interactions within 
the chains, the magnetic relaxation is expected to present 
an unusual behavior. In particular, slow dynamics have 
been investigated in single chain magnets^ as well as 
in ordered spin chains^. Recently, unexpected resonant 
effects, attributed to quantum tunneling, have been ob- 
served in the three-dimensional (3D) ordered spin chain 
[(CH 3 )3NH]CoCl3-2H 2 (CoTAC) 7 .. 

The study of avalanches in such systems can provide 
additional information about the reversal processes. A 
pertinent question concerns the starting flip of a popu- 
lation of spins at the origin of the avalanche. A quan- 
tum origin was proven in the early results on the Mn^ac 
single molecule magnet (SMM) 8 which appears as an 
emblematic example of intrinsic avalanches in molecular 
compounds. The avalanches were shown to occur around 
the resonance field and were thus attributed to resonant 
quantum tunneling of the magnetization which increases 
the spin nipping rate^. 

Another important question concerns the propaga- 
tion of the avalanche and thus the associated time 
scale required for a macroscopic magnetization reversal 
to develop. For example, recent avalanches studies in 
Mni2ac have shown that the magnetization reversal is 
not uniform inside the sample 9 and that the avalanches 



propagate at a constant velocity, requiring a threshold 
energy-^. In most cases, the propagation of the avalanche 
can be attributed to thermally assisted phenomena. At 
low temperature, it is often difficult for the heat released 
by the spin flipping process to be dissipated in the sam- 
ple and to be absorbed by the external thermal bath (due 
to poor thermal coupling of the sample to the environ- 
ment and/or low thermal conductivity). In this scenario, 
the local heating due to the flipping of a small group or 
cluster of spins is sufficient to heat the neighboring spins 
to a temperature where thermal activation is efficient, 
thus enabling them to flip, which in turn will heat their 
neighbors, and so on. 

To know whether or not the latter mechanism governs 
the avalanche, two experimental tests can be done. One 
is simply to vary the sample size and the coupling to 
the heat bath. Indeed, the use of quite small samples 
allows better thermal homogeneity and avoids excessive 
overheating during a magnetization reversal process, thus 
suppressing the avalanche. For example, when a small 
crystal of Mn^ac was directly immersed in liquid 4 Hc- 
3 He inside the mixing chamber of the dilution refriger- 
ator, the sudden avalanches were replaced by anomalies 
in the hysteresis loop (called steps) , which are the signa- 
ture of the relaxation of the magnetization by quantum 
tunnclingii. Another test is to vary the field sweeping 
rate. During a field ramp, a certain amount of small 
stochastic spin flips may occur, which are due to de- 
faults or impurities in the sample and may be consid- 
ered extrinsic. However, if enough are present and if the 
field is ramped too fast, then the combined heat released 
from all of these processes during a short time can trig- 
ger a thermal avalanche. The field at which this kind 
of avalanche occurs will depend on the temperature and 
ramping rate in a very complicated way. On the other 
hand, by ramping more slowly, the heat will have time to 
dissipate and the avalanches will disappear. This was the 
case at low temperature for CoTAC 7 -, where, for example, 
avalanches occurred at field values as low as 100 Oe when 
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the field was ramped at 75 Oe.s^ 1 . By simply ramping 
at a slower rate of approximately 1 Oe.s -1 , avalanches 
were suppressed, and we were able to measure relaxation 
in fields up to 2000 Oe and, in particular, explore the 
resonant tunneling near 1000 Oc. 

In the present paper, we are interested in the avalanche 
phenomena in a metallo-organic system, a Mn (Il)-free 
radical ferrimagnetic chain. Contrary to the above sys- 
tems, here, we show that the avalanches are intrinsic 
both in their origin and their propagation. Indeed, the 
avalanches are found to be sample independent at low 
temperature and to occur at a well-defined internal field. 
They do not disappear even by using very slow field 
sweeping rates, down to 0.04 Oe.s^ 1 , and with samples 
as small as 10 -4 mm 3 . 

In Sec. [Ill we briefly describe the samples, review 
some pertinent experimental features, and summarize the 
structural and magnetic properties from previous inves- 
tigations on this compoun d 12 ' 13 . In Sec. IIIIl we present 
our magnetic measurements at very low temperature (be- 
low 800 mK), focusing on the avalanches detected in the 
hysteresis loops. Since these avalanches appear to be in- 
trinsic to the system, we propose in Sec. IIVI that the 
observed behavior is related to domain-wall propagation. 
We present a simple model that allows us to describe the 
main experimental features. 



II. SAMPLES AND MAGNETIC 
CHARACTERISTICS 

The studied compound is the Mn(hfac) 2 (i?)-3MLNN 
(formula C^TtbsF^Mn^Og) 1 ^ comprised of chains of 
Mn and nitronyl-nitroxide free radicals (NITR), called 
here Mn-NITR. We obtained our main results for three 
samples of different shapes and weig hts (see Table H}. 
Samples 1 and 2 are needle shaped and were mea- 
sured using low temperature SQUID magnetometers 
equipped with a miniature dilution refrigerator. The 
set-up can measure absolute values of the magnetization 
by the extraction method 1 ^. Sample 3, about a thou- 
sand times smaller, was measured with a micro-SQUID 
magnetometer—. All the measurements we present here 
were performed along the easy axis, the b axis. In this 
direction, the values of the demagnetizing factor N are 
quoted in Table [J They were deduced from magnetiza- 
tion measurements, sample shape and, in case of samples 
1 and 2, ac susceptibility. 



Sample Mass Dimensions (mm) N (cgs) 
Sample 1 1.68 mg 4.2x0.95x0.95 ^1 
Sample 2 0.818 mg 4.7x 0.6 x0.5 ~0.6 
Sample 3 <1 fj,g 0.1x0.05x0.05 ~2.5 



TABLE I: Summary of measured samples. 

The magnetic structure of the compound is described 
in a separate paper-i 3 -. It originates directly from the 



crystallographic structure 12 . The zig-zag chains extend 
along the 6-axis and are composed of alternating chiral 
NITR free radicals ((i?)-3MLNN=(i?)-Methyl[3-(4,4,5,5- 
tctramethyl-4,5-dihydro-l//-imidazolyl-l-oxy-3-oxide) 
phenoxy]-2-propionate)) carrying a spin snitr = 1/2 
and Mn(hfac)2 units where the Mn(II) ion carries a 
spin ^Mn = 5/2. The spins preferentially point along 
the chain axis. This point, which is unusual in Mn(II) 
and radical based compounds, was adressed in Ref. 
Il3lwhere it was shown that the anisotropy cannot be 
explained by dipolar interactions: Single- ion anisotropy 
or antisymmetric exchange has to be invoked to account 
for the anisotropy^. The interaction between the Mn 
and radical spins along the chains is antiferromagnetic, 
very strong, Jintra ~ 500 K, and much larger than the 
interchain coupling J int er (Jintra/ ' J inter ~ 10 4 ). The 
magnetic susceptibility follows a purely ID ferrimagnetic- 
like behavior at high temperature, crossing over to a ID 
ferromagnetic behavior below 90 K involving effective 
spins Stot = Smh — snitr = 2. An analysis of the 
susceptibility as a function of temperature in the range 
of 20 K< T <90 K allows us to estimate an equivalent 
ferromagnetic interaction J = 65 K between the effective 
spins Stor^-- 

Below 20 K, interchain coupling effects become impor- 
tant, and give rise to a 3D long-range magnetic order be- 
low T c ~ 3 K. Magnetization and ac susceptibility mea- 
surements above T c suggested the existence of a ferro- 
magnetic transition at T c . On the other hand, magneti- 
zation measurements performed below T c failed to show 
a spontaneous magnetization. Neutron diffraction mea- 
surements removed this contradiction. In zero field, the 
order is actually antiferromagnetic and is composed of 
alternating planes of ferromagnetic aligned chains. How- 
ever, the order is tenuous, and a small field is sufficient 
to push the system into the ferromagnetic state (i.e., ~ 
150 Oe at 1.6 K)^ 3 -. This behavior has been interpreted 
in terms of a subtle competition between the interchain 
correlations which are of the order of a few milliKclvins: 
long-range dipolar interactions which favor an antifer- 
romagnetic order vs. weak short-range super-exchange 
interactions, responsible for ferromagnetic correlations 
(their presence was also inferred from measurements un- 
der pressure 1 ^): As the temperature is decreased, corre- 
lations along the chains grow, and the dipole field be- 
comes proportionally stronger. At T c , it overwhelms the 
weak super-exchange interactions and 3D antiferromag- 
netic order is established. 



III. RESULTS: DYNAMIC PROPERTIES 

A. Hysteresis loops and magnetic avalanches 

We focus here on the dynamic properties in the 3D 
magnetically ordered phase below 800 mK, at temper- 
atures where hysteresis loops occur. From measure- 
ments of the magnetization performed along the chain 
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FIG. 1: M vs. H: Decreasing part of the hysteresis loop for 
sample 1 at several temperatures between 140 and 400 mK 
at an equivalent field sweeping rate of 0.4 Oe.s -1 (steps of 1 
Oe every 2.5 s). The sample was saturated in a 2000 Oe field. 
Avalanches appear below 300 mK. 
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FIG. 2: Coercive field H c o as a function of temperature for 
sample 3 at a field sweeping rate of 350 Oe.s - . Empty circles: 
No avalanches occur. Full circles: H c o is equal to the first 
avalanche field. The line is a guide to the eyes. The insets 
show representative hysteresis loops in both regimes. 



axis (which is the easy magnetization direction), we per- 
ceive two distinct temperature regimes, readily apparent 
in the hysteresis behavior shown in Figs. Q] and [2] In the 
high temperature regime (above 300 mK), the magneti- 
zation M(H) decreases smoothly as the magnitude of the 
applied (negative) field increases. Below 300 mK, a new 
physics develops and one observes magnetic avalanches. 
In Fig. [TJ we show the decreasing part of the hystere- 
sis loops for a few illustrative temperatures for sample 1. 
The data are obtained after the application of a positive 
2000 Oe field sufficient to saturate the magnetization 18 . 
In the following, we define the coercive field H co as the 



field at which M(H) crosses the M = axis, whether 
avalanches are detected or not. The smooth hysteresis 
curves observed at high temperature give way to abrupt 
avalanches below 300 mK. Note that one avalanche is 
sufficient to totally reverse the magnetization at T > 200 
mK but several avalanches are needed for the same pur- 
pose below 200 mK (See Fig. [T]). Nevertheless, in the 
latter case, the amplitude of the first avalanche was al- 
ways sufficiently large so as to cross the M = axis, so 
conveniently H co is also the threshold field for occurrence 
of the first avalanche. The coercive field depends on T, 
but not monotonously, in the studied temperature range. 
See the examples of sample 1 in Fig. [T]and of sample 3 in 
Fig. O When decreasing the temperature, one observes 
that H co first increases down to 300 mK (regime with 
no avalanches), subsequently decreases down to 200 mK 
and finally becomes temperature independent below 200 
mK. Here, we stress that the hysteresis features are really 
intrinsic, and sample independent. In fact, although the 
sample shapes and masses varied, and different appara- 
tuses and procedures were used for the measurements, we 
always found a maximum H co occurring at 300 mK. Also, 
below 300 mK, we always observed magnetic avalanches 
for all samples, regardless of the value of the field sweep- 
ing rate (See below). 

Once an avalanche takes place, it typically lasts for ap- 
proximately 15 ms for the 1 mg Mn-NITR samples. This 
leads to a 10-15 m.s -1 propagation velocity, comparable 
to Mni2ac avalanche velocity2ii£. We could also detect 
a faint temperature pulse during the avalanche of a few 
millikclvins (the thermometer was approximately 10 cm 
away). This is, by comparison, much smaller (at least 
ten times) than what we previously observed in similar 
size Mni2ac or CoTAC samples (taking into account the 
different avalanche field values). We will see below that 
the characteristic internal field at which avalanches take 
place is -ffi=170 Oe. Therefore, we can estimate the 
heat released for an effective spin flipping in this field: 
AT = gno^BASHx/CLT = 1.9 K (in which the low tem- 
perature specific heat Clt ~ 4 x 10 6 erg.mol -1 .K -1 was 
measured by the quasi-adiabatic method down to 350 
mK— ). Each effective spin is surrounded by six first 
neighbor spins located on adjacent chains. As a chain 
flips, we imagine that the heat radiates away from the 
chain like spokes on a wheel. In first approximation, if 
all of this heat is absorbed by the nearest chains, it would 
result in a local heating that corresponds to a tempera- 
ture rise of only w 300 mK per spin. When compared 
to the characteristic energy scales of the system (see Sec. 
IIII B|) , this seems to indicate that local heating is not the 
only mechanism that contributes to the avalanche pro- 
cess. 

A crucial result concerns the value of the internal field 
Hi aval for the first avalanche. Hi ava i is obtained from 
the applied field H corrected for demagnetization effects 
in a mean field approach: Hi ava i = H — NM , where N 
is the demagnetization factor noted in Table U and M 
is the value of the magnetization just before the occur- 
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FIG. 3: First avalanche field Hi aval, corrected for demag- 
netizing field, as a function of field sweeping rate v, for two 
different temperatures, 100 and 200 mK. Squares, respectively 
circles and crosses correspond to sample 1, respectively sam- 
ple 2 and 3. Lines indicate the field Hi = 170 Oe. 



rence of the avalanche. At 100 mK, and when the sam- 
ple is in the saturated state, we found the same value of 
Hi = Hi avai{Mb = M sat ) = 170 Oe for all samples what- 
ever their shape or volume is (varying by a factor 10 3 ) 
and independent of the field sweeping rate w, that we var- 
ied by 4 orders of magnitude from 0.04 to 700 Oe.s -1 (See 
bottom of Fig. [3]). In addition, the value of Hi is very 
reproducible: When repeating hysteresis loops about 200 
times for samples 1 and 2, we found the first avalanche 
field to be constant to within 2 %. Furthermore, as seen 
in Fig. [TJ below 200 mK, the reversal occurs in several 
avalanches. It is worth noting that, when corrected for 
demagnetization effects, the successive avalanches occur 
at the same internal field of 170 Oe (but with a larger 
distribution than the first avalanche field). These prop- 
erties lead us to conclude that this avalanche field is an 
intrinsic property of Mn-radical chains. 

At 200 mK, Hi aV ai is still constant, but only for sweep- 
ing rates faster than 0.5 Oe.s -1 . For field sweeping rates 
slower than 0.5 Oe.s -1 , Hi ava i is no longer constant (Sec 
top Fig. [3|). In fact, the experimental conditions have 
changed when the field ramping is very slow at this higher 
temperature. This is because the magnetization has time 
to relax during the sweep due to thermally activated pro- 
cesses and therefore the magnetization M just before the 
avalanche differs substantially from the saturation value. 
This suggests that the value of the magnetization just 
before the avalanche M plays a crucial role in the deter- 
mination Of Ht aV al ■ 

To clarify this point, we performed a series of M(H) 
measurements at 100 mK, starting from a non-saturated 
state. For these measurements, we field cooled the sam- 
ple in various fields from above 2 K (where there is no 
hysteresis) down to 100 mK in order to freeze the sam- 
ple in different given magnetization configurations, with 
different amounts of initial magnetization. The result- 
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FIG. 4: The first avalanche field Hi aval vs. the initial magne- 
tization Mf, just before an avalanche for sample 2, at 100 mK 
and for field sweeping rates of 9.6 Oe.s -1 (full circles) and 
0.12 Oe.s" 1 (empty circles). The inset shows M(H) curves 
obtained at 0.12 Oe.s -1 for several Mi, values. 



ing avalanche field Hi a vai at 100 mK as a function of 
the magnetization just before an avalanche M is shown 
in Fig. Q] for two field sweeping rates. As can be seen, 
Hi aval increases with decreasing initial magnetization, 
and seems to saturate at about 300 Oe. Note, in par- 
ticular, that the avalanche field starting from the zero 
field cooled (ZFC) magnetization H2 = Hi ava \(M b = 
0) « 270 ± 20 Oe2£ is larger than the avalanche field 
of Hi = 170 Oe obtained from the saturated state (See 
Fig. 01 . One can see why a constant avalanche field is 
observed at low temperature when starting from satura- 
tion: As the temperature is decreased, relaxation during 
the hysteresis loop becomes increasingly slow. Below 300 
mK, the starting value for M remains very close to sat- 
uration, and according to Fig. 31 H t aV ai will occur at 
approximately 170 Oe. 



B. Estimation of the characteristic energy barrier 

The rapid variation of H co with temperature between 
300 and 900 mK suggests an origin in terms of thermally 
activated processes. The latter can be well characterized 
from the study of the dependence of the magnetization 
on the measuring time t (or frequency /) and tempera- 
ture T. Thus, we studied the dissipative part x" [T, f) of 
the ac susceptibility (H ac \\b) for samples cooled in zero 
static field (i.e. in the antiferromagnetic phase). \" (T, /) 
exhibits a maximum at a temperature Tm that depends 
on the measuring frequency and obeys an Arrhenius law: 
t = To exp(i?o/fcslM), t = l/2_7r/ over six decades in 
frequency (See Fig. [5] and Ref. [l2j ). The energy barrier 
Eq is sample independent: For samples 1, 2 (but also 
for two other samples not presented here), we found the 
same value of Eq/Ub = 16 K even though the samples 
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FIG. 5: (Color online) lnr vs. 1/T for four different samples, 
from ac susceptibility measurements with H a c ~ 0.14 Oe and 
1.1 mHz< / < 511 Hz (r = l/27r/). For each sample, the line 
represents the fit to an Arrhenius law: r = to exp(_Bo/fcsT), 
with the energy barrier Eo/ks « 16 K and 7 x 1CP 11 s < to < 
5 x 10 -10 s. The inset shows the dissipative part x" of the 
ac susceptibility as a function of temperature at frequencies 
between 0.57 and 111 Hz for a 1.59 mg single crystal. 
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FIG. 6: Distribution of energy barriers, P(E). Squares (cir- 
cles): distribution obtained from relaxation vs. time t at 
several temperatures on sample 1 (sample 2), Triangles: dis- 
tribution obtained from relaxation vs. temperature T of sam- 
ple 2. The inset shows the magnetization M as a function of 
E/k B = Tln(t/ro) with to = 2.8 x 10~ 10 s from relaxation 
curves as a function of time of sample 2. All these relaxation 
curves were measured in zero field. 



were synthesized in different batches over a two year pe- 
riod and had various thermal cycles. Only tq is found 
to vary slightly between samples, ranging from 7 x 10~ n 
to 5 x 10~ 10 s , and may be explained by different char- 
acteristic chain lengths from sample to sample^. Cole- 
Cole plots of isothermal x"(/) vs - x'(f) are non-circular 
which implies a distribution of energy barriers P(E) of 
which Eq is the mean value. An analysis of our x" (u>) 
and data using models for the distribution of en- 

ergy barriers^, provides a rather narrow width of P(E), 
of about AEi/2/kB = 3 K at half maximum^. 

Below 800 mK, x"(T) vanishes. However, it is still 
possible to investigate Eq and P(E) from 800 mK down 
to 300 mK by studying the relaxation of the saturated 
isothermal remanent magnetization Mr in zero field. 
Note that during this relaxation process, the system 
passes from the field induced ferromagnetic phase to the 
antiferromagnetic ground state. We did not find a sim- 
ple exponential relaxation characteristic of the presence 
of a single barrier height. Instead, we found that the re- 
laxation of Mr depends on a single variable, E c (T,t) = 
fcsTln(t/To) and is characteristic of a distribution of bar- 
rier heights crossed by thermal activation. In this case, 
P(E) can be found by adjusting the Mu(T,t) data to 
the formula P(E = E c ) = dMa/dE^. Recording the 
evolution of the magnetization at different constant tem- 
peratures over 10 h following the removal of the field, 
we could superpose the Mr(T, t) data on a single curve 
in the scaling diagram Mr vs. Tln(t/ro) when setting 
t = 2.8 x 10~ 10 s (See the inset of Fig. EJ. From the fit, 
we deduce a P(E) centered on the energy E /kB =15.5 
K with /S.E l / 2 /k B =3 K (See Fig. O, thus in very good 
agreement to that deduced from our ac data in the anti- 



ferromagnetic phase above 800 mK. Thus, using different 
experimental methods (in different temperature ranges), 
we found nearly the same energy barrier Eq and distri- 
bution P(E), independent of the sample used. 

Finally, we note that we also measured, after saturat- 
ing the sample, the relaxation of the magnetization in a 
negative applied field. Below 300 mK, and for field values 
less than Hi, we observed two regimes. At short times 
(beyond our experimental resolution), the magnetization 
shows a small decrease (few percent). At long times (up 
to 24 h), the relaxation is very slow, with a logarithmic 
tail, and at 100 mK less than a fraction of a percent of 
the magnetization relaxes. 



IV. DISCUSSION 
A. Origin of the energy barrier 

In the following, we propose that the avalanche dy- 
namics at low temperature are governed by the same en- 
ergy barrier observed at higher temperature, which we ar- 
gue arises from the energy required to nucleate a domain 
wall along a chain. However, the avalanche dynamics are 
found to be modified by the state of the magnetization 
Mb just before the avalanche. We suggest that it is a con- 
sequence of the competition between ferromagnetic and 
antiferromagnetic orders in Mn-NITR, and propose a sce- 
nario that allows us to describe the observed behavior by 
considering two situations: the avalanches from the sat- 
urated ferromagnetic state and those from the ZFC an- 
tiferromagnetic state. Finally, we discuss the questions 
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that emerge from this scenario and note that the local 
heating and the subsequent diffusion of the heat front, 
although important, are not the driving forces behind 
the propagation of the avalanche. 

We contend that the energy barrier E /kB ~ 16 K ob- 
served in measurements above 300 mK is the energy re- 
quired to nucleate a domain wall in the system. Although 
Mn-NITR orders three dimensionally and there are very 
strong ID magnetic interactions within the chains, the 
interactions between spins on neighboring chains (which 
ultimately result in the 3D order), are very weak, only 
of the order of a few millikelvins (See Sec. |TTJ) . Because 
the interchain interactions are so weak compared to the 
anisotropy along the chain, there can be no domain wall 
between chains. That is to say, for a given chain, the 
neighboring chains will either be parallel or anti-parallel, 
and thus for directions perpendicular to the chains, there 
is no "width" that one usually associates with a domain 
wall. 

On the other hand, as already evidenced in other spin 
chain compound s 7 ' 247 , along the chains the competition 
between the exchange and anisotropy energies results in 
the energetically favorable conditions for the nucleation 
of 0D domain- walls. To estimate the energy of such a 
domain- wall, we will, for simplicity, consider in this pa- 
per that the chains are made of effective spins S to t — 2. 
This is justified since all our experimental data below 
20 K 12 i 13 ' 19 are consistent with an effective spin approx- 
imation. Futhermore, as the NITR spin is delocalized, 
the spin distribution between the Mn and the NITR is 
not knowni 3 -, making it very difficult for a more detailed 
description of the chain. Then, the cost in exchange en- 
ergy at low temperature needed to flip an effective spin 
Stot — 2 on a chain edge would be JSj ot = 260 K, a very 
large energy. The cost to flip a spin inside a chain would 
be twice as much. However, the nucleation of a 0D do- 
main wall along a chain will have a much smaller activa- 
tion energy. In a first approximation, we consider a 180° 
domain-wall with energy 2 ^: 

JSl^ nK 

2n 2 v ; 

where n is the number of spins Stot in the domain- 
wall, J is the exchange constant between them, and K is 
the anisotropy constant, supposing a uniaxial anisotropy. 
Eq. [T] describes a domain- wall in which each spin is ro- 
tated by an angle of n/n radians with respect to its neigh- 
bors. The first term is the cost in exchange energy for 
such a domain- wall, which decreases as 1/n and thus fa- 
vors a more spread-out wall. The second term is the 
increase in anisotropy energy for such a domain-wall and 
is proportional to n. 

Magnetization isotherms and perpendicular suscepti- 
bility allow us to estimate the anisotropy constan t 13 ' 19 : 
K « 2 x 10 4 erg. cm -3 , that is to say 0.13 K/spin Stot- 
The domain-wall energy F (See Eq. [1} is minimized for 



n = y — — = 140 effective spins in the chain. This 

rather large wall is a consequence of the weak anisotropy 
in this compound. The energy required to nucleate the 
wall is F/ks = 18 K. This value is very closed to the ex- 
perimental energy barrier Eo/ks — 16 K notwithstand- 
ing the simplicity of the model. Local defaults can modify 
this nucleation energy, which can explain the observed 
barrier distribution. Once a wall is nucleated, it will 
sweep along the chain reversing many more spins. 

We can use this rough estimate of the number of spins 
in a domain wall in order to compare the experimen- 
tal energy barrier with the Zeeman energy acquired by 
the wall just before avalanche takes place. Consider for 
example, the case when the sample has first been sat- 
urated in high field, and is thus in the ferromagnetic 
state. Avalanches always occur at H\ = 170 Oe when 
Mb = M sat . Assuming that the same number of spins 
(ri = 140) is involved in the nucleation of the domain- 
wall, the Zeeman energy supplied by the field to these 
spins is Eu 1 /kB = gfisnStotHi — 6.4 K (with g = 2 
from EPR measurements 27 ). In a similar fashion, when 
the sample has been zero field cooled in the antiferro- 
magnetic state with M = 0, avalanches take place when 
H% w 270 Oe. In this case, the Zeeman energy acting to 
flip the spins and create the wall is approximately 10.1 
K. The above estimates of the energy scales, although 
systematically less than the experimental values, are nev- 
ertheless of the same order of magnitude. 

Because the ferromagnetic state in zero field has a 
higher energy than the antiferromagnetic ground state 
in Mn-NITR, we might expect that it is easier to nucle- 
ate a domain wall in the ferromagnetic state. This may 
explain in part the dependence of the Zeeman energy on 
the initial magnetization. Indeed, it has been showni 3 - 
that dipole interactions are responsible for the 3D anti- 
ferromagnetic phase transition at 3 K. So it seems rea- 
sonable to estimate the difference between the two states 
from magneto-static energy considerations. At H = 0, 
the magnetic energy density in the ferromagnetic state is 
-1/2 M sat Hi) = 740 erg. cm -3 (where Ho is the demag- 
netizing field with N = 1). This corresponds to an in- 
crease of approximately 4.5 mK/spin 2 above the ground 
state. For a domain wall of 140 spins, this implies a 
difference in energy of about 0.65 K for the wall, signif- 
icantly smaller than the observed difference of 3.7 K in 
the Zeeman energy for the walls. Intriguingly, however, 
this is of the same order of magnitude as the difference 
between the energy barriers measured in the antiferro- 
magnetic state by the ac susceptibility and the slightly 
smaller barrier deduced from relaxation out of the anti- 
ferromagnetic state. 



B. Phenomenological model 

We suggest a simple model to describe qualitatively the 
observed behavior with the aid of Fig. [7] which schemat- 



ically shows the energy barrier landscape for two differ- 
ent initial states. In Fig. [7^, the sample has been first 
saturated in a high field and is in the ferromagnetic- 
up state at H = 0. The ferromagnetic-up state and 
ferromagnetic-down state are degenerate at H — 0, but 
separated by an energy barrier. The antiferromagnetic 
ground state lies slightly below these. At low tempera- 
ture, thermal activation is not enough to overcome the 
energy barrier and the system remains in the excited fer- 
romagnetic state. Fig. [TJd shows the effects of applying a 
magnetic field. In a negative field, the Zeeman energy of 
the ferromagnetic states is shifted with respect to each 
other: The up state increases in energy, the down state 
decreases by an equal amount, while the antiferromag- 
netic state remains the same. If the temperature is very 
low during the field ramp, relaxation is very slow and Mb 
remains constant and close to the saturation value M sat . 
At the critical field Hi , the energy barrier of the up state 
will be reduced to near zero, and multiple nucleations of 
domain walls will occur, creating an avalanche. The sys- 
tem will slide down into the lowest state, which, due, to 
the applied field, is now the ferromagnetic down state. 

The situation when the sample has been zero field 
cooled in the antiferromagnetic Mb = ground state 
is shown in Fig. [7J:. In this larger field must 

be applied in order to reduce the energy barrier and in- 
duce nucleations of domain walls, as shown in Fig. (7J1. 
The system then avalanches into the ferromagnetic-down 
state. Note that, in this simple model, it has been as- 
sumed that nucleation of domain walls occur when the 
energy barrier is reduced to nearly zero. Another pos- 
sibility could be that the energy barrier remains finite, 
and nucleation of the domain- walls takes place by quan- 
tum tunneling through the barrier. This would, in effect, 
short circuit the energy barrier and explain why our esti- 
mates for the Zeeman energy of the walls are lower than 
the experimental energy barrier. 

It is interesting to look at the multiple avalanches at 
low temperature in light of our simple model. As seen 
in Sec. IIII A[ when corrected for demagnetization ef- 
fects, the successive avalanches are found to occur at 
nearly the same internal field H\ « 170 Oe. This re- 
sult implies that the sample remains in the ferromag- 
netic state during the avalanche process, but is broken 
up into a few large ferromagnetic domains of opposing 
polarity. This is in agreement with the scenario shown 
in Fig. [TJo. When avalanching out of the M sat state, 
the sample goes from the ferromagnetic-up state to the 
ferromagnetic-down state directly; i.e., it does not pass 
through the antiferromagnetic state. This also implies 
that the internal field is critical to the advancement of 
the ferromagnetic domain front. During the avalanche, 
the internal field is reduced below its critical value as the 
net magnetization changes. If the avalanche stops, then 
at the interface of two oppositely polarized domains, the 
field will be a minimum. Then, the external field has 
to be increased until the critical internal field is reached 
again, so that a new avalanche is induced. 



H=0 Mb=M sa t 



En 



j i H=H1 
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II 



II 



H=H2>H1 
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(ci) 



II 



FIG. 7: Schematic diagrams of the energy barrier landscape, 
a) Mb = M 3a t and H — 0: The system has been field cooled in 
the ferromagnetic-up state. The nucleation energy barrier E n 
has to be crossed to go from the ferromagnetic-up state to the 
antiferromagnetic ground state, b) Mb = M 3a t and H — Hi: 
Because of the Zeeman energy, the ferromagnetic-down state 
has the lowest energy, and Hi is enough to suppress the bar- 
rier between the ferromagnetic-up and down states. The sys- 
tem avalanches into the ferromagnetic-down state, c) Mb = 
and H — 0. The system has been zero-field cooled in the anti- 
ferromagnetic state which is separated from the others by the 
energy barrier E' n . d) M b = and H = Hi > Hi: The field 
is enough to suppress E' n between the antiferromagnetic state 
and the ferromagnetic-down state. The system avalanches 
into the ferromagnetic-down state. 



This simple model does not address the apparent col- 
lective phenomena of the avalanche. This aspect is quite 
different from those previously observed in Mn^ac SMM 
or the CoTAC chain. As the latter is also 3D ordered^, it 
is interesting to compare: As mentioned in the Introduc- 
tion, avalanches could be suppressed in CoTAC samples 
by simply ramping the field at a slower rate. At very low 
temperature, the relaxation of the magnetization from 
+M sat to — M sat could be measured. The cause of the 
relaxation was shown to be due to quantum nucleation of 
domain walls (made up of only ten spins or so) that took 
place at a resonant field of approximately 1000 Oe. Each 
nucleation of a domain wall by tunneling through an en- 
ergy barrier could be treated as an independent event. 

This is very different from the present situation for Mn- 
NITR: At low temperature, avalanches always took place 
regardless of the sample size or slow ramping speed. At 
100 mK, the relaxation of the magnetization was always 
extremely slow, r > 10 6 s, right up to Hi ava i. Of course, 
when an avalanche took place, it was rapid, t rs 10~ 2 
s. After an avalanche, there was no further relaxation, 
which excludes the hypothesis of a resonance effect as ob- 
served in CoTAC. Thus, there is an all or nothing aspect 
to the dynamics at low temperature of Mn-NITR. Other 
important differences between the two systems include- 
the following: the spins in the Mn-NITR chains come 
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from two different sources (the Mn(II) ion and from the 
delocalized electron of the free radical NITR) that are 
aligned along the chain direction, the ID character is 
more than an order of magnitude stronger than in Co- 
TAC, and the anisotropy is more than an order of mag- 
nitude weaker. 

These attributes result in an estimated domain wall 
length that is more than an order of magnitude longer 
in Mn-NITR in comparison with CoTAC. This leads us 
to speculate that the importance of the long-range dipo- 
lar interchain interactions in Mn-NITR may provide the 
conditions for the collective nucleation of domain walls. 
Within the long and spread-out domain wall, a large 
number of spins will be perpendicular to the chain axis. 
This will produce a sizable transverse field on the neigh- 
boring chains, which, in turn, should favor the induce- 
ment of a domain wall in the neighbors. 

That local heating operates in this system but does 
not dominate the avalanche process can be seen in the 
following illustration. Below 200 mK, avalanches always 
occurred in a few large, distinct steps, as can be seen in 
Fig. Q] for T = 140 mK. The magnetic energy released 
when large blocks of spins flip, 2MHi ava i is converted 
to heat during the avalanche. However, at low tempera- 
ture, this heat is not enough to raise the temperature high 
enough for thermal activation to be efficient and, thus, 
not enough to sustain the avalanche and finish the sample 
off. On the other hand, above 200 mK, avalanches were 
always complete, that is to say, from +M sat to —M sat in 
one step. Thus, the same magnetic energy of 2MHi ava i 
generated during the avalanche, along with the additional 
thermal energy due to the higher starting temperature, 
was enough to overwhelm the sample and cause a com- 
plete magnetization reversal. 



V. CONCLUSION 



In conclusion, we have shown that intrinsic avalanches 
occur in Mn-NITR, which is unusual in macroscopic sam- 
ples. Below 200 mK, they are characterized by a constant 
reversal field in the hysteresis loops, which does not de- 
pend on the sample size or batch, the temperature and 
the field sweeping rate. We propose that this avalanche 
field is related to the energy barrier measured in ac sus- 
ceptibility and relaxation measurements and corresponds 
to the energy needed to nucleate domain- walls along the 
chains. We suggest that, by applying a field, we tilt the 
energy field landscape to the point that the barrier height 
goes to zero, and due to the purity of the sample and the 
sharp energy distribution, many domain- wall nucleations 
occur, resulting in an avalanche. We also speculate that 
this is a collective phenomena. 
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